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Abstract 

Equations are derived which describe a propagation of strong shocks 
in the interstellar matter , without any demands for a symmetry, in a 
thin layer approximation (2.5 dimensions). Using these equations permits 
to calculate a propagation of shock waves from nonsymmetric supernovae 
explosions in the medium with arbitrary density distribuion, formation of 
superbubbles in galaxies. 



1 Thin layer approximation 

A thin shell approximation for discription of strong shocks is based on two 
simplifications. First, it is assumed that all swept-up intercloud gas accumulates 
into the thin shell just behind the shock front and moves with the velocity u. 
Second, the pressure distribution inside the cavity P{r, t) is taken to be uniform. 
The equations of mass and momentum conservation may be expressed as follows 
(Cherny, 1957) : 

M = Mo + 47r r pir^dr, 4^-R'(Pin - P) + Mg, (1) 

Jo dt 

where M is the mass of the shell, Mq is the ejected mass, R is the shock radius 
and u is the gas velocity behind the shock; p(r) and P are the density and 
the pressure of the ambient gas, g is the external gravitational field. For the 
adiabatic blastwave without gravity 

"f = " ^'f ^^o = E.„4m„=, ,2) 

where Eq — const is the energy of the explosion, Eth = ^(^^zrrj PinR^ is the 
thermal energy of the blastwave, and 7 is the adiabatic power. Equations ( |^)- 
( H) have a simple solution for the homogeneous case if the swept-up mass is 
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much greater than the ejected one R = 
Bisnovatyi-Kogan and BUnnikov, 1982) 



Po 



1/5 



^2/5, where (Chernyi, 1957, 



75(7-l)(7+l)' 



2)^0 = 



75(7 - 1) 



1/5 



(3) 



167r(37- 1) 

Comparison with an exact self-similar (SS) solution (Sedov, 1946) shows, that 
the case 2) gives much better precision, ^0^^ = 1.033, 1.014 1.036 for (SS), (1), 
(2) cases respectively at 7 = 1.4; and similarily = 1.15, 1.12, 1.15 for 
7 = 5/3. 



2 Three-dimensional shocks 

Introduce following Bisnovatyi-Kogan and Silich (1991), Silich (1992) (see also 
Palous, 1990, Bisnovatyi-Kogan and Sihch, 1995) Cartesian coordinate system 
{x,y, z). Let m is the mass, r is the radius- vector, u is the velocity of the par- 
ticular Lagrangian element of the shock, p{x, y, z) = pof(x, y, z) is the ambient 
gas density, n is the normal to the shock front unit vector, g is the acceleration 
of the external gravitational field, V is the velocity field of the undisturbed gas 
flow, S is the surface area of the Lagrangian element, m = ctS, cr is the surface 
density, AP = Pi„ — P is the pressure difference between the hot interior and 
warm (cold) external gas. The pressure Pin = (7 — l)Eth/^ of the hot tenuous 
gas within the cavity is a function of the bubble thermal energy Eth and vol- 
ume fl. To describe shock expansion we must introduce element of the surface 
area and define volume of any closed three-dimensional region. It is well known 
from the differential geometry that any surface may be specified parametrically, 
with the Cartesian coordinates at any point on the surface being function of 
two parameters, Ai and A2: x = x{Xi, A2), y = y{Xi, A2), z = z{Xi, A2). Then 
the element of the surface area may be defined by the expression (Budak and 
Fomin, 1965) 



= 5(Ai,A2)(iAidA2, 
d{y, z) 



5(Ai,A2) = 



9(Ai,A2 



d{z, x) 



a(Ai,A2) 



+ 



d{x,y) 



1/2 



(4) 
(5) 



where d{xi, Xj) / d{\i, X2) = det{d{xi, Xj)/d{Xi, X2)) are the appropriate Jaco- 
bians. Similarly could be obtained expressions for Q and n 

If parameters Ai and A2 are considered as the Lagrangian coordinates of the 
shock front, then the equations for 3-D shock propagation may be written for 
the mass /x = cr S'(Ai,A2) per unit of Lagrangian square, in a compact form, 
convenient for numerical integration 
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d/i 

dux AP d{y, z) - 



■PX + 5x, (7) 

PX + 9v, (8) 

+ (9) 

h- (10) 



di yu 9(AiA2) /i 

dUy AP d{z,x) Uy — Vy 

At fi d{XiX2) jJL 

du^ AP d{x, y) - 



dt ^ d{XiX2) ^ 
dx dy dz 



Here function x is defined as follows: 



/ d{y,z) d{z,x) d{x,y) 

Approximating the shock front by a number N of Lagrangian elements one get 
a system of 7N differential equations of the mass and momentum conservation. 
This set of equations is coupled with the equation for the gas pressure within 
the cavity, and the equation of energy for Etot — Eth + Pk + Eg, consisting from 
the thermal energy of the hot bubble interior, kinetic and gravitational energies 
of the shell. The kinetic and gravitational energies of the shell are determined 
by corresponding surface integrals. Variations of the total energy Ptot of the 
remnant throughout the adiabatic stage of evolution are defined by the energy 
input rate L{t), kinetic and thermal energies of the swept- up interstellar gas 
with temperature T{x, y, z): 



Etot =En+ [ [Lit) + ' """^ / ' """^ + 3fcr/?7)dAidA2]dt, (12) 

where Eq is the initially deposited energy, rj is the mean mass per particle. On 
the radiative phase of expansion the gas behind the shock front cools so quickly 
that does not add to the total energy of the remnant. Rarefied hot gas inside 
the cavity expands adiabatically and accelerates the surrounding dense shell. 
Time-derivative of the thermal energy of the remnant is defined then by the 



equation, which is used instead of ( 12) 
dEth 



dt 



L{t)- / Pi„u„5(Ai,A2)dAidA2, (13) 
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where u„ is the velocity component normal to the shock front. Calculations 
of supershell formation in the plane-stratified differentially rotating Galactic 
disk have been performed by Silich et al. (1994) and are represented in the 
Fig. ^ An hourglass remnant with a noticeable degree of deformation by the 
Galactic shear have appeared. Formation of elongated superbubbles due to 
differential galactic rotation gives a possibility to determite unambiguously a 
direction of the galactic rotation (Mashchenko and Silich, 1995). The means of 
this determination is shown in the Fig. ||. 

Figure 1: Galactic superbubble morphology for different locations of the parent 
OB-association relative to the galactic plane. Left: the OB-association is at the 
midplane of the Galaxy. Right: the OB-association is 50pc above the Galactic 
plane. 

Figure 2: Scheme of the HI holes orientation in the plane of view. 
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